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Abstract 

We consider the integrable open XX quantum spin chain with nondiagonal bound- 
ary terms. We derive an exact inversion identity, using which we obtain the eigenvalues 
of the transfer matrix and the Bethe Ansatz equations. For generic values of the bound- 
ary parameters, the Bethe Ansatz solution is formulated in terms of Jacobian elliptic 
functions. 



1 Introduction 



We consider in this article the open XX quantum spin chain with nondiagonal boundary 
terms, defined by the Hamiltonian 



N-l 

^ = 2{E(«+l + <^n+l) 
n=l 



+ ^cothe-ai^ + — ^a^-^cothe+a^ + — ^a^j, (1.1) 

where (7^,(7^, are the standard Pauh matrices, and ^± , k± are arbitrary boundary pa- 
rameters. This model is known to be integrable |]T], 0, |^. It has been investigated Q 
using a fermionization technique p, |^, suitably adapted to accommodate boundary terms 
P). However, this model has until now resisted a direct Bethe Ansatz solution due to 
the absence of a simple reference (pseudovacuum) state. Q Such a solution is desirable for 
a number of reasons. First, the open XX chain is a special case of integrable open XXZ 
and open XYZ chains, which should also admit Bethe Ansatz solutions but which cannot 
be solved by fermionization. Second, Bethe Ansatz solutions are particularly well-suited 
for investigating physical properties, such as ground state, low-lying excitations, scattering 
matrices, etc. In particular, the Bethe Ansatz approach avoids the projection mechanism 
which can be implemented only in special cases. (See 0, Q and references therein for discus- 
sions of interesting physical applications of the open XX spin chain.) Finally, the Sklyanin 
transfer matrix for the open XX spin chain is closely related to the Yang matrix [|ll|, [T^, |13[ 
for a large class of integrable M = 2 supersymmetric quantum field theories with boundary 



14, 15]. Diagonalization of this matrix is a key step in formulating the thermodynamic 



Bethe Ansatz equations for these Af = 2 supersymmetric models. 



In this paper we derive an exact inversion identity for the model (|1.1| ), using which 
we obtain the eigenvalues of the transfer matrix and the Bethe Ansatz equations. (This 
approach does not completely circumvent the problem of not having a reference state, since 
the eigenvectors are not determined.) We obtain the inversion identity using the open-chain 
fusion formula |jl6[ together with the remarkable fact that, for the open XX spin chain/ 6- 
vertex free-Fermion model, the fused transfer matrix is proportional to the identity matrix. 
A similar strategy has recently been used [|T^ to solve the open 8-vertex free-Fermion 
model, which corresponds to the case of A/" = 1 supersymmetry. These techniques are 
generalizations of those which have been developed for closed spin chains [|l^, |l^, |2T|| . 

Even though the transfer matrix is constructed entirely from hyperbolic functions, we find 

^For the special case of diagonal boundary terms (i.e., k± = 0), a simple pseudovacuum state does exist, 
and a Bethe Ansatz solution is known 
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that the Bethe Ansatz solution is formulated in terms of Jacobian elliptic functions for generic 
values of the boundary parameters. (In contrast, for the XYZ chain |jl8|, such functions 
appear already in the transfer matrix.) For special values of the boundary parameters, the 
elliptic functions degenerate into ordinary hyperbolic or trigonometric functions. 

The outline of this paper is as follows. In Section 2, we review the construction of the 
Sklyanin transfer matrix for the open XX chain, and derive some of its important properties. 
We derive in Section 3 the inversion identity, which we then use in Section 4 to determine the 
Bethe Ansatz solution. In Section 5, we investigate some special cases in which the solution 
can be expressed in terms of ordinary hyperbolic functions. In particular, we verify that 
our solution is similar to the known one |jlO|, |l[ for the case of diagonal boundary terms. In 
Section 6, we conclude with a brief discussion of some possible directions for future work. 



2 Transfer matrix 



The object of central importance in the construction of integrable quantum spin chains is 
the one-parameter family of commuting matrices called the transfer matrix. The transfer 
matrix for an open chain is made from two basic building blocks, called R (bulk) and K 
(boundary) matrices. 

An R matrix is a solution of the Yang-Baxter equation 

Ruiu - v) RM RM = R23iv) RM R^iu - v) . (2.1) 

(See, e.g., [^, |2^.) The XX spin chain is a special case of the XXZ spin chain, corre- 
sponding (in the 
the 4x4 matrix 



spending (in the notation of [||]) to the anisotropy value = y. The R matrix is therefore 



R{u) 



/ a 








o\ 





b 


c 








c 


b 





V 








a) 



(2.2) 



with matrix elements 



a = sinh(ti -|- 



b = sinh u , c = sinh 



iir 



which satisfy the free-Fermion condition 

2 I 1,2 2 

a + = c 



(2.3) 



(2.4) 
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This R matrix has the symmetry properties 

Ri2{u) = Vi2Ri2iu)Vu = RuiuY''' , (2.5) 

where V12 is the permutation matrix and t denotes transpose. Moreover, it satisfies the 
unitarity relation 

Ruiu) Rui-u) = C{u)I , C{u) = - cosh^ u , (2.6) 
and the crossing relation 

Ru{u) = ViRu{-u-pyWi, (2.7) 

with 

p = --, V = a\ (2.8) 
Finally, it has the periodicity property 

Ruiu + irt) = -alRi2{u)al = -alRr2{u)al . (2.9) 

The matrix K~ (u) is a solution of the boundary Yang-Baxter equation 

Ru{u - v) K^{u) R2i{u + v) K^{v) = K^{v) R^iu + v) K^{u) R2i{u - v) . (2.10) 

We consider here the following 2x2 matrix 0, ^ 

^_^^^^/sinh(e- + t.) /._sinh2« \ ^^.ll) 
\ K_smh.2u sinh((^_ — ?i) / 

which evidently depends on two boundary parameters , We set the matrix K^{u) to 
be K^{—u — p) with (^_ , k_) replaced by ((^+ , k+); i.e., 

^^^^^^/zcosh(e^-«) .^sinh2. \ ^2.12) 
\ K+sinh2M — z cosh(,^+ + m) / 

We shall often use an alternative set of boundary parameters (77^= which is related 

to the set (^^p , k^:) by 

i 1 

cos ?7=p cosh = sinh f=r , cos^ 7?=r + cosh^ = 1 ^ ^ . (2.13) 

The K matrices have the periodicity property 

K^{u + m) = -a'K^{u)a\ (2.14) 



The transfer matrix t{u) for an open chain of spins is given by |l| 

t{u) = tioK^iu) Toiu) K,{u) foiu) , (2.15) 

where tro denotes trace over the "auxihary space" 0, and Tq{u), Tq{X) are so-called mon- 
odromy matrices 

To(ti) = Ron{u) ■ ■ ■ Roiiu) , fo{u) = Rioiu) ■ ■ ■ Rm{u) . (2.16) 

Indeed, Sklyanin has shown that t{u) constitutes a one-parameter commutative family of 
matrices 

[t{u),t{v)] = 0. (2.17) 

Typically, the Hamiltonian is proportional to the first derivative of the transfer matrix t'{0) 
|l[]. However, this quantity is trivial for the XX model, due to the fact tri^+(0) = 0. In 
order to obtain the Hamiltonian, we must go to the second derivative We find 

t"(0) 



n 



4(_l)^+i^sinh tiK+'iO) 



iV-1 

2sinh^_"^ ' tri^+'(0) ' ^ ^ ^ 

n=l ^ ' 

where Hn,n+i = Vn,n+iR'n n+i{^) ^ ^^^1 wc have made use of the facts i^^(O) = sinh^_I and 
troK'(^(0) Hno = trfC^"(0) = 0. By explicitly evaluating ( p.l8| ), we obtain the Hamiltonian 
(1.1). Notice that the Hamiltonian is Hermitian if ^± are imaginary and k± are real. The 



corresponding energy eigenvalues E are given by 

E = ^ (219) 

8(-l)^+isinhe-sinh^+ ' ^ ' 

where A(m) are eigenvalues of the transfer matrix. 

The transfer matrix has the periodicity property 

t{u + m) = t{u) , (2.20) 

as follows from ( |2.9D , ( |2.14| ). Moreover, the transfer matrix has crossing symmetry 

t{-u - '-^) = t{u) , (2.21) 

which can be proved using a generalization of the methods developed in the appendices of 
26| . Finally, we note that the transfer matrix has the asymptotic behavior (for k± ^ 0) 

g«{4+2Af) 

t{u) ^ K^K+i^ 1+ . . . for u^oo. (2.22) 
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As is customary, we usually suppress the "quantum-space" subscripts 1 , . . . ,N. 
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3 Inversion identity 



Our main objective is to determine the eigenvalues A(tt) of the open-chain transfer matrix 
( p.l5| ), from which the energy eigenvalues (p.l9|) immediately follow. We shall accomplish 
this using an exact inversion identity, which we first derive. A similar approach was used in 
T9| for closed chains. This approach is based on the concept of fusion pO|, 



The derivation of the inversion identity for the open XX spin chain/ 6- vertex free-Fermion 
model closely parallels the one for the 8- vertex free-Fermion model considered in |13|. For 
brevity, we shall often refer to these two models as the M = 2 and M = 1 models, respectively. 
The principal tool which we use to derive the inversion identity is the open-chain fusion 
formula obtained in 11161. We shall henceforth refer to this reference as I. 



The matrix Ru^u) at u = — p = y is proportional to the one- dimensional projector 



12 2 



(^1 



12) 



P. 



12 



(3.1) 



12 

/ \ 
110 
110 

V / 

As explained in I, from the corresponding degeneration of the (boundary) Yang-Baxter equa- 
tion, one can derive identities which allow one to prove that fused (boundary) matrices satisfy 
generalized (boundary) Yang-Baxter equations. 

The fused R matrix is given by (I 2.13) 

R<12>3{U) = P+ RM R23{U + p) P+ , (3.2) 

where Pj^ = I — Pi2- An important observation is that the fused R matrix can be brought 
to the following upper triangular form by a similarity transformation 



Xi2 R 



<12>3 



U) X 



-1 
12 



f Sl * 

ta' 

V 







* \ 



0/ 



(3.3) 



where 

s = cosh^ u , 
and the 4x4 matrix X is given by 



t = i cosh u sinh u , 



(3.4) 



X 



/ 1 \ 

10 1 

1 

\ 1 1 / 



(3.5) 
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It follows that the fused monodromy matrices (I 4.7), (I 5.4), (I 5.5) 

T<i2>iu) = R<i2>n{u) ■ ■ ■ R<:i2>i{u) , 
T<i2>{u + p) = R<i2>i{u) ■ ■ ■ R<i2>n{u) , 
also become triangular by the same transformation. 



Xl2 T^12>{u) 



= Xi2 T<i2>(n + p) Xi2\ 

where F = nii<- 

The corresponding fused K matrices are given by (I 3.5), (I 3.9) 
^<i2>(w) = Pi\Ki{u) Ri2{2u + p) K^{u + p) P,\, 
Kti2>{n) = {Pt2 Ktiur Ri2i-2u - 3p) K^iu + p)*^P+}*" 

since M = V^V = I. 



(3.6) 







* 


* 


*\ 







t^F 


((-1)^- l)t^F 













(-t)^F 





v 











0/ 



(3.7) 



(3.8) 



Unlike the A/" = 1 case [13], the similarity transformation does not bring also the fused 
K matrices to upper triangular form. Q Nevertheless, the transformed fused K matrices are 
"almost" triangular 



X,2 K^i2>iu) X{,' 



* \ 







777-4 '^s 



4 

V / 



(3.9) 



where 



rrin 



mi 



ml 



±- sinh 2u [cosh 2u — cosh 2^zp + sinh^ 27i] , 



— - sinh 2u [sinh 2 (w ± ^^) — cosh u sinh^ u\ 



- sinh 2u [sinh 2 (w =F ^^) + cosh^ u sinh u] 



in^ cosh 2u sinh 2u , 



±- sinh Au sinh 2^=p . 



(3.10) 



■^It is not possible to simultaneously triangularize both i?<i2>3(u) and -K'<i2>(u), since their commutator 
is not nilpotent; i.e., [-R<i2>3(u) ,^<i2>(w)]" 0. A monograph on the general problem of simultaneous 
triangularization has recently been published p7| . 
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The fused transfer matrix t{u) is given by (I 4.5), (I 4.6) 

i{u) = tri2 K^i2>{u) T<i2>{u) i^<i2>(^i) ^<i2>(m + p) ■ 



(3.11) 



With the help of the results (|3^ ), (|3^ ),( ^30| ), we now obtain the remarkable result that the 
fused transfer matrix is proportional to the identity matrix, 



t{u) = A{u)I, 



(3.12) 



where 



AW 



4 ^5 



+ (1 — (—l)^)(m3 777-4 + ^3 ~ m^m^ — m2 ml + 27774 ^4 ) 

A similar result also holds for the M = 1 case ||l3l . 
The fusion formula is given by (I 4.17), (I 5.1) 
1 



(3.13) 



t{u) t{u + p) 



C(277 + 2p) 



i{u) + A {K+{u)]/\ {K-{u)]5{T{u)} 6 \ T( 



u 



(3.14) 



where the transfer matrix t{u) is given by (|2.15| ) (see also (I 4.1), (I 4.2)), and the quantum 
determinants ^ are given by (I 4.15), (I 5.3), (I 5.7) 



6{T{u)} = 6[T{u)\=au + pf, 
A{K-{u)} = tTu{PuK^{u) Ri2{2u + p) K^{u + p)ViV2} , 
A {K+{u)} = tri2 {Pfa V, K+{u + p) R,2{-2u - 3p) K+{u)} 



(3.15) 



It follows from ( ^.12| )-( P.15| ) that the transfer matrix obeys an exact inversion identity 



t{u) t{u--) = f{u)I. 



(3.16) 



where the function f{u) is given by 

f{u) = tanh^ 2u [gi{u) cosh"^^ u + g2{u) sinh^^ u + gz{u) sinh^^ u cosh^^ , (3.17) 

with 

Qi {u) = - (cosh 2u — cosh 2,^_ + sinh^ 277) (cosh 277 — cosh 2^+ + 2/€^ sinh^ 27i) 
= 1Qk\h^_ cosh(77 + 77^_) cosh(t7 — iriJ) cosh(7t + irjj^) cosh(77 — 77^+) 

X COsh(M + 7?_) COsh(M — d-) COSh(77 + 7?+) COsh(77 — 7?+) , 
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92{u) 



2k\k^_ [I + (-1)^ + sinh^ 2u] sinh^ 2u + (-1)^ \{k\ cosh + cosh 2^+) sinh^ 2w 



+ -(cosh2,^_ cosh2(^_|_ sinh^ 2w — sinh 2^_ sinh2^+ cosh^ 2m) 



2 1 t\j 



sinh^ 2u cosh^ 2u + (— 1)^ ^ sin 2?7_ sinh 2'(9_ sin 2?]+ sinh 2'(9+ cosh^ 2u 



+ cos 2?7_ cosh 2-d_ cos 277_|_ cosh 2-19 _|_ sinh 2w 



(3.18) 



The inversion identity ( P.16| )-( p.l8| ) is the first main result of our paper. We have checked 
it numerically up to = 3. 



4 Eigenvalues and Bethe Ansatz Equations 

Having obtained the inversion identity, we now use it to determine the eigenvalues of the 
transfer matrix. The commutativity relation (|2.171) implies that the transfer matrix has 
eigenstates |A) which are independent of m, 



t{u)\A) = A{u)\A) 



(4.1) 



where A(m) are the corresponding eigenvalues. Acting on |A) with the inversion identity, we 
obtain the corresponding identity for the eigenvalues 



ITX , 



A{u) A{u--) = f{u). 



(4.2) 



Similarly, it follows from ( |2.20| ) and ( |2.21| ) that the eigenvalues have the periodicity and 
crossing properties 



A{u + in) = A{u) , A{-u ) = A{u) . 

2 



(4.3) 



Finally, ( p.22[ ) implies the asymptotic behavior 

^ g«(4+2Af) 



N' 



+ ... 



for u ^ oo . 



2l+2iV 

We shall assume that the eigenvalues have the form 

N 

A{u) = psinh 2u JJ^ sinh(w — Uj) cosh('U + uj) 

j=0 



(4.4) 



(4.5) 
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where Uj and p are (m- independent) parameters which are to be determined. Indeed, this 
expression satisfies the periodicity and crossing properties ( [4 .31) , and it has the correct asymp- 
totic behavior ( [4 .41) provided that we set 

p = i^4:K_K+ . (4.6) 



We now substitute the Ansatz ([4.5|) into the inversion identity (f4.2| ), and obtain 



^ 1 

l)^p2 sinh^ 2m JJ - sinh 2(m - Uj) sinh 2{u + uj) = f{u) . 



(4.7) 



j=0 



Recalhng the exphcit expression ( |3.171 ),( pT8D for f{u), we verify that both sides of the 
equation have the same asymptotic behavior ~ g«(8+4Af) ^ qq_ Since the LHS has zeros 
±Uj, these must be zeros of f{u). That is, Q 



gi{uj) cosh"^^ Uj + g2{uj) sinh"^^ Uj + gsiuj) sinh^^ Uj cosh^^ uj = . 
Dividing by cosh^^Wj, we obtain 

g2{uj) tanh^^ Uj + g^iuj) tanh^^ Uj + gi{uj) = . 
Regarding ( [4.9| ) as the quadratic equation 

g2x'^ + g3X + gi = Q 



(4.9) 



(4.10) 



in the variable x = tanh^^ Uj, we conclude that the parameters Uj satisfy the Bethe Ansatz 
equations 



tanh^^ Uj 



h{u 



1/ 



92[Uj 



(4.11) 



where the function h{u) is defined by 



h{u) 



-gsiu) ± A/^^(wp^^%I(ti)~^ 



u 



(4.12) 



The square root in ( [4.12| ) can be eliminated by making an appropriate change of variables. 



Indeed, with the help of ( 3. IS ), one can show that 

h(u) = K^K^_< — sinh^ 2?/ cosh^ 2-u — (— 1)^ (71 cosh^ 2^ + 72 sinh^ 2^) 



dz sinh 2u cosh 2u \j a sinh^ 2u + {3 1 , 



(4.13) 



*The three functions gi{u) ,i = 1,2,3, are even functions of u. Hence, if Uj is a root of f{u), then so is 
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where 

a = - (cos477„ + cos 477+ + cosh4'(9_ + cosh4'!9+) + (— l)^[cos2(?7_ + 77+) cosh2('i9_ — 
+ cos 2{r]^ — 77+) cosh 2('!9_ + , 

2 1 

P = [sin277_sinh279_ + (-I)^sm2?7+sinh27?+] - -[sm2(r/_ + 77+) sinh2(i?_ - 

+ sin2(?7_ — ?7+) sinh2('(9_ + 79+)]^ , 

71 = sin 2?7_ sinh 2'(9_ sin 27/+ sinh 2^9+ , 

72 = cos 2?7_ cosh 2'i9_ cos 277+ cosh 2i?+ . (4-14) 
Let us change from the spectral parameter u to the new spectral parameter v defined by 

sinh 2m = z sn 277 , cosh 27i = cn 2f , (4-15) 
where the modulus k of the Jacobian elliptic functions is given by 

k' = ^. (4.16) 
With the help of the identities (see, e.g., p9[] ) 



cn^;2 + sn^z = 1, dn^ z + A;^ sn^ ^ = 1 , (4.17) 

one can see that the function h{u) can be reexpressed as 

h{u) = K^K^I sn^27;cn^27; - (-1)^ (71 cn^ 2t; - 72 sn^ 2t;) ±i sn 27; cn 27; dn 2f a/^I (4.18) 

Hence, the Bethe Ansatz solution (|]|), (U), ( pl| ) can be reformulated as 

1 ^ 
A(u) = {--)^-\+K_sn2vYl{sn2v - sn2vj) , (4.19) 



j=0 

where the parameters Vj satisfy 



cn 27;^ — 1 \ ^ sn^ 27;^- cn^ 27;j — (—1)^ (71 cn^ 27;^ — 72 sn^ 27;^) ± 7 sn 27;^- cn 27;^ dn 2vj\f^ 



cn 277^ + 1 / (cn 27;^ — cn 27r/_) (cn 2vj — cn 279_) (cn 2vj — cn 2277+) (cn 27;^- — cn 27?+) 

(4.20) 

This Bethe Ansatz solution is the second main result of our paper. This result passes 
several tests. Indeed, for = 0, 1, the eigenvalues agree with those obtained by direct 
diagonalization of the transfer matrix. Moreover, as discussed in the following section, our 
solution is similar to the known one 0, Q for the case of diagonal boundary terms. 
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5 Special cases 



For generic values of the boundary parameters, the Bethe Ansatz solution presented in the 
previous section is formulated in terms of Jacobian elliptic functions. However, for modulus 
k = or k = 1, these elliptic functions degenerate into ordinary trigonometric or hyperbolic 
functions. Equivalently, the argument of the square root in ([4.13|) then becomes a perfect 
square, and so the square root effectively disappears. We now briefly consider some of these 
special cases. 



5.1 Diagonal case 



In the limit k± — >■ 0, the K matrices ( |2.11| ), ( |2.12| ) become diagonal, and therefore so do 



the boundary terms in the Hamiltonian ( |1 . 1|) . The transfer matrix t{u) now commutes with 
the operator F = YliLi^i^ hence, both operators can be simultaneously diagonalized. 
Denoting the corresponding eigenstates by lA'-^-'), we have 

t(tx)|A(±)) = A(±)(ii)|AW), 

F|A(±)) = ±|AW). (5.1) 

The transfer matrix now has the asymptotic behavior 

^(^) ~ P 21+2N ^ + • ■ • u-^oo, (5.2) 



where 



. z cosh((^_|_ — (^_) for N= even (5 3) 

z sinh((^+ — ,^_) for N= odd 



The eigenvalues are given by 



N 

A^^\u) = ipsinh 2tt JJ^ sinh(tt — Uj) cosh(M + Wj) , (5.4) 



which is similar to ([4.5|) except with one less root. 



In terms of the boundary parameters ri± (|2.13|) , the limit k± ^ corresponds to 



TT 1 

77± = - - , e"'^* = ^ oo . (5.5) 

In this limit the function h{u) becomes equal (for = even) to 

h{u) = — sinh(M =]= cosh(ii =]= ^_) sinh(M ± ^+) cosh(M ± ^+) . (5.6) 
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Moreover, the Bethe Ansatz equations ( [4.11|) become 



cosh(wj ± (^_) cosh{uj =F ^+) 



(5.7) 



These resuhs are similar to those obtained previously |1T0| , |l| for the XXZ chain. The two 
choices of signs correspond to the two possible pseudovacua - either all spins up or all spins 
down. 

5.2 Nondiagonal cases 

Within the space of boundary parameters, there are various submanifolds, such as 

V--V+± - = ^^^^^ ' (5-8) 

for which a = k = 0. 

As a simple example, let us consider the particular case N = odd, rj^ = rj^ = r], = 
= 'd, for which also j3 = 0. The function h{u) becomes equal to 

h{u) = —K^K^ sinh 2(m — irj) sinh 2(m + irj) sinh 2(m — -i?) sinh 2{u + i?) . (5.9) 

We then obtain the Bethe Ansatz equations 

tanh^^ Uj = — coth(Mj + irj) coth('Uj — it]) coth(wj + i!)) coth(Mj — {)) . (5.10) 

6 Discussion 

This work raises a number of interesting questions, some of which we list below: 

We have seen that the doubly-periodic functions in the Bethe Ansatz solution degenerate 
into singly-periodic functions for special values of the boundary parameters. In ||^, important 
simplifications are also found to occur for special values of the boundary parameters. It is 
likely that these two observations are related. 

As we have emphasized, an exact inversion identity for the XX (or M = 2) case is 
made possible by the key fact that the fused transfer matrix is proportional to the identity 
matrix. A similar result also holds for the M = 1 case [|1^. It would be interesting to better 



understand the relation of this phenomenon to the free-Fermion condition. 

The model ( |1 . 1| ) is not the most general integrable open XX chain. Indeed, the most 
general solution 0, 0] of the XXZ boundary Yang-Baxter equation has the off-diagonal terms 
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K^"* sinh 2u and k,^2^ sinh 2u, while here we have restricted to the special case k^2^ = k^'' = k±. 
(See Eqs. ( p.llD , ( p.l2D .) However, we do not expect that the more general case will lead to 
new significant complications. In particular, we expect that the same approach can be used 
to derive an exact inversion identity and to obtain the corresponding Bethe Ansatz solution. 



Since the model (|1.1| ) has various boundary parameters, its phase diagram is likely to 
have a rich structure. Our exact Bethe Ansatz solution should provide a means of exploring 
these phases. Moreover, as mentioned in the Introduction, this solution opens the way to 
formulating the thermodjTiamic Bethe Ansatz equations for integrable M = 2 supersym- 
metric quantum field theories with boundary |n, 15|, as was done for the case of A/" = 1 
supersymmetry in |^ . Finally, with the insight gained from the XX chain, it might now be 
possible to finally solve the open XXZ chain with nondiagonal boundary terms. 

We hope to report on some of these problems in future publications. 
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